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ABSTRACT: 

In  an  investigation  of  the  cooling  system  for  the  ORNL  Research 
Reactor,  calculations  have  been  made  to  relate  the  pressure  loss  in 
the  assemblies  with  the  flow  rate.  In  addition  an  expression  was 
developed  and  evaluated  to  give  the  maximum  temperature  of  the  surface 
in  contact  with  the  water  as  a  function  of  flow  rate,  inlet  water 
temperature,  and  maximum  heat  flux.  After  this  maximum  heat  flux  has 
been  determined  and  a  maximum  permissible  surface  temperature  selected, 
the  corresponding  minimum  flow  rate  and  its  associated  pressure  loss 
can  be  found  for  any  inlet  water  temperature  by  the  methods  illustrated 
in  this  report.  Examples  of  such  results  are  shown  in  Table  V. 

As  a  by-product  of  the  pressure  loss  calculations,  certain 
recommendations  have  been  made  pursuant  to  the  reduction  of  pressure 
loss  in  the  assemblies.  This  recommendation  is  that  the  inlet  and 
exit  flow  areas  be  made  rectangular  in  cross  section,  thereby  in¬ 
creasing  the  flow  area  and  eliminating  any  sudden  change  from  cir¬ 
cular  to  rectangular  cross  section. 
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HEAT  TRANSFER  AND  PRESSURE  LOSS  IN  PROPOSED  ORR  FUEL  ASSEMBLIES 

INTRODUCTION 

In  order  to  design  the  most  efficient  cooling  system  for  the  ORNL 
Research  Reactor,  it  is  necessary  to  know  the  relationships  between  flow 
rate  through  the  fuel  assemblies,  the  pressure  loss,  and  the  maximum 
surface  temperature  at  various  operating  powers.  Since  no  flow  tests 
have  been  made  for  the  proposed  ORR  fuel  assemblies,  the  pressure  loss 
as  a  function  of  flow  rate  must  be  calculated  from  empirical  flow  re¬ 
lationships.  The  maximum  surface  temperature  can  be  calculated  from 
a  simple  heat  balance  if  the  assumption  is  made  that  no»heat  is  con¬ 
ducted  linearly  along  the  fuel  plate. 

CALCULATION  OF  PRESSURE  LOSS? 

Description  of  the  Assembly- 

Ess  entially,  the  assembly  is  similar  to  a  standard  MTR  fuel 
assembly  with  the  exception  that  round  adapters  are  used  on  both  ends 
so  that  the  element  may  be  inverted  for  more  complete  burn  up.  The 
assembly  is  described  in  ORNL  Drawing  D-15898  with  the  exception  that 
two  modifications  have  been  made.  These  changes  are  (1)  it  was  as¬ 
sumed  that  entrance  and  exit  straight  sections  shown  as  2.000  inches  I,  D. 
could  be  opened  to  2.2^0  inches  I.  D.  and  (2)  the  recent  changes  de¬ 
scribed  in  ORNL  CF  53-12-119^  in  the  MTR  fuel  and  side  plates  are  in¬ 
corporated  in  this  design.  These  changes  include  making  the  fuel  plates 
50  mils  thick  which  include  a  20  mil  plate  of  uranium-aluminum  alloy 
clad  on  both  sides  with  l5  mil  aluminum  sheets.  The  side  plates  are 
reduced  in  thickness  from  3/l6  inch  to  1/8  inch. 


*  Cunningham,  J.  E. ,  ^Information  Paper  on  MTR  Fuel  Elements,”  ORNL  CF 
53-12-119,  December  21,  1953. 
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In  reducing  the  fuel  plate  thickness  from  the  original  design 
of  60  mils,  the  overall  cross  sectional  dimensions  of  2.996  inches  and 
3.069  inches  are  maintained.  If  the  original  number  of  plates,  18,  is 
maintained,  this  increases  the  gap  width  between  plates  to  0.1276  inches 
from  the  original  0.117  inches.  If  an  additional  plate  is  added,  the 
gap  width  becomes  0.1177  inches,  practically  the  original  dimension.. 

Since  there  was  some  doubt  as  to  the  advantages  of  adding  this 
nineteenth  fuel  plate,  calculatioifs  have  been  made  using  both  configu- 
rations  in  the  pressure  loss  calculations. 

Empirical  Relationships  Used  to  Determine  Pressure  Loss- 

The  losses  due  to  sudden  contractions  were  calculated  according 
to  Equation  (31)  in  Badger  and  McCabe.*  Changing  this  relationship  from 
head  loss  to  pressure  loss  in  pounds  per  square  inch,  it  becomes 

K  /o  (1) 

c  / 

^  P  "  2  dull)  gc 

The  constant,  K,  is  evaluated  according  to  Figure  1$  of  this  reference.** 
c 

The  significance  of  the  remaining  terms  is  found  in  Appendix  I  of  this 
report.  The  losses  due  to  a  gradual  contraction  are  neglected.*** 

The  losses  due  to  sudden  expansions  were  calculated  from 
Equation  (30)  in  Badger  and  McCabe<***«  which  is  identical  with  Equation 
(102)  in  Vennard.'tt^HttHt 


"Badger,  Ifl/.  L.  and  McCabe,  W.  L.,  Elements  of. 
Second  Edition,  McGraw-Hill  Book  Co.,  Inc., 


Chemical  Engineering, 
New  York,  193^^  V*  hP* 


^l^ard*,  J.*K.,  Elementary  Fluid  Mechanics,  Second  Edition,  John  Wiley 
and  Sons,  Inc.,  New  York,  19U^*  P*  17U. 

Badger,  W.  L. ,  and  McCabe,  W.  L. ,  0£.  £it .  >  P*  37. 

Vennard,  J.  K. ,  o£.  cit . ,  p.  17U. 
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Again  changing  to  pressure  loss  in  pounds  per  square  inch,  the  equation 


(Ujj  -  u  )2^ 

AP  "  _ 

2(1UU)  Sc 

In  the  case  of  a  more  gradual  expansion,  the  formula  is  modified  according 
to  Equation  (103 )  in  Vennard,^* * * **** *****“ 

“  2(l!iU)  gc 

The  constant,  K^,  was  evaluated  according  to  Figure  8l|  of  this  reference, -JHt 

The  losses  due  to  friction  in  straight  lengths  of  constant  cross 

sectional  area  were  calculated  according  to  the  Fanning  friction  equation 

shown  in  McAdains.*«*  The  integrated  form  of  Equation  (8a)  of  this  text, 

converted  to  pressure  loss  in  pounds  per  square  inch,  is 

2  f  u^/’l 
“  mrg'c  Be ' ' 

The  Fanning  friction  factor,  f,  is  found  hy  computing  the  Reynolds  number, 
then  applying  the  relationship  stressed  graphically  in  Figure  $1  of 
McAdams. A  line  drawn  proportionately  one-fourth  the  distance  from 
line  A  to  G  on  the  plot  was  found  to  best  relate  the  friction  loss  in  the 
flow  between  the  fuel  plates. 

In  applying  the  friction  equation  and  in  determining  the 
Reynolds  number  for  flow  in  channels  other  than  circular,  and  eouivalent 
hydraulic  diameter  was  calculated  according  to  the  suggestion  in  Badger  and 
McCabe, The  equivalent  diameter  is  given  as  four  times  the  cross- 
sectional  area  of  the  channel  divided  by  the  wetted  perimeter  of  the  channel. 


*  ^id,  p.  175. 

•tHf  Ibid,  p,  176. 

McAdams,  W,  H, ,  Heat  Transmission,  Second  Edition,  McGraw-Hill  Book  Co,, 
Inc.,  New  York,  (19U2},  p.  119. 

****  Ibid,  p,  118. 

*****  Badger,  W.L. ,  and  McCabe,  W,  L. ,  0£.  clt. ,  p.  i^||. 
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Method  of  Calculation- 

In  the  calculation  of  the  pressure  loss  through  the  fuel  assembly, 
the  velocity  between  the  fuel  plates  was  assigned  the  value,  V*  By  calcu¬ 
lating  the  relative  flow  areas  throughout  the  assembly,  the  velocities  at 
other  points  could  be  determined  in  terms  of  this  symbol*  The  pressure  loss 
over  each  section  was  then  calculated,  also  in  terms  of  this,  symbol.  In 
order  that  the  pressure  loss  calculated  might  coirpare  with  the  pressure  loss 
which  might  be  determined  in  an  experimental  apparatus,  the  effect  of  changes 
in  velocity  head  was  calculated.  The  net  change  in  velocity  head  over  the 
entire  assembly  was  zero  since  it  was  assumed  that  velocities  were  the  same 
above  and  below  the  assembly  support  grids. 

Calculations  were  made,  as  was  stated  previously,  both  for  18- 
and  19-fuel  plate  assemblies.  These  calculations  were  made  assuming  that  the 
average  water  temperature  was  lOO^F,  For  the  19— fuel  plate  assembly,  additional 
calculations  were  made  assuming  an  average  water  temperature  of  160^,  The 
viscosity  of  water  at  these  teir^eratures  was  determined  from  the  nomograph. 
Appendix  II,  in  Badger  and  McCabe.*  The  density  of  water  at  these  tem¬ 
peratures  was  determined  as  the  reciprocal  of  the  specific  volume  given  in 
Keenan  and  Keys,** 

An  example  of  these  calculations  is  shown  in  Appendix  II  of  this 

report. 

*  Ibid.,  p.  632. 

**  Keenan,  J.  H, ,  and  Keys,  F.  G.,  Thermodynamic  Properties  of  Steam, 

First  Edition,  John  Wiley  and  Sons,  Lric.,  iew  lork,  (19l;8J,  p.  29-30, 
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Validity  of  these  Helationships- 

In  order  to  check  the  validity  of  the  application  of  the  formulae 
to  the  situation,  a  calculation  of  the  pressure  loss  in  a  standard  MTR  fuel 
assembly  was  made  and  compared  to  experimentally  determined  results.  The 
pressure  loss  was  calculated  for  an  average  linear  water  velocity  of  thirty 
feet  per  second  between  the  fuel  plates  at  an  average  temperature  of  80^F, 
These  were  compared  with  experimental  data  given  in  ORNL-CF-50-6-102.* 

The  results  are  shown  in  Table  I, 

TABLE  I 

COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL  PRESSURE  LOSS 
IN  MTR  FUEL  ASSEMBLIES  FOR  A  WATER  VELOCITY  OF  30  FT/SEC 


Calculated 

(Psi) 

Experimental 

(psi) 

Upper  End  Box 

lli.l 

13 

Fuel  Plate  Section 

17.1 

18 

Lower  End  Box 

9.0 

6 

Total 

itO.2 

37 

^  Kasper,  S.,  "Water  Flow  in  Materials  Testing  Reactor,"  ORNL  CF  50-^-102, 
June  1$,  1950>-.p.  6. 
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Results  of  Pressure  Loss  Calculations- 

The  relationship  giving  the  pressure  loss  through  the  assembly  as  a 
function  of  velocity  through  the  fuel  plate  is  given  in  the  following  express- 


ions. 


In  the  18-plate  assembly  at  100°F, 

Upper  end  box  •••••••« 

..  Ap 

Fuel  plates  . 

.•  Ap 

Lower  end  box  . 

..  Ap 

Total  Assembly  .. 

..  ^P 

In  the  19-plate  assembly  at 

0 

0 

0 

H 

Upper  end  box  . . 

.. 

Fuel  plates  . .  •  •  • 

..  AP 

Lower  end  box  . . 

AP 

Total  Assembly  ... 

Zip 

In  the  19-plate  assembly  at  160^F, 


Upper  end  box  . .  AP 

Fuel  plates  Ap 

Lower  end  box  . . . .  •  •  AP 


Total  Assembly  AP 


»  0.0092^  +  0.00109 

-  0.00387  +  0.0272 

-  0.01012  +  0.00195 

o  1  79 

-  0.0232U  T  +  0.0302  V  * 

2  1.79 

=  0.00857  V  +  O.OOlOU  V 

-  0.0039U  +  0.0291 

-  0.00960  +  0.00186 

-  0.0221  +  0.0320 

-  o.ooShU  +  0.00098 

-  0.00388  +  0.0276 

1  79 

«  0.009UU  V  +  0.00175  V  ' 

-  0.0218  +  0.0303 


1 
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The  velocity,  V,  is  given  in  feet  per  second,  and  the  pressure  loss  is  in 
pounds  per  square  inch.  Table  II  gives  values  for  the  pressure  loss  at 
various  velocities.  Figure  1  shows  this  relationship  graphically.  Exam¬ 
ination  of  this  plot  indicates  that  the  variation  in  total  pressure  loss 
between  these  three  situations  does  not  differ  significantly. 


TABLE  II 

PRESSURE  LOSS  IN  ORR  FUEL  ASSEMBLIES 
with  18-  and  19-  plates 
as  a  function  of 
Velocity  between  the  Plates 


Velocity 

V 

(ft/sec) 

Pressure  Loss  in 
19-plate  Assemblies 
at  lOO^F  (psi) 

Pressure  Loss  in 
18-plate  Assemblies 
at  100°F  (psi) 

Pressure  Loss  in 
19-plate  Assemblies 
at  16oOF  (psi) 

1 

0.05 

0.05 

0.05 

2 

0.20 

0.20 

0.19 

1* 

0.77 

0.77 

o.iu 

8 

2.7i| 

2.7ii 

2.65 

10 

U.19 

U.19 

i;.05 

15 

9.06 

9.06 

8.77 

20 

15.69 

15.76 

15.10 

25 

2U*03 

21t.l6 

23.28 

30 

3U.0 

3U.3 

3U.0 

35 

1i5.7 

U6.0 
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CALCULATION  OF  KAXIMJM  SURFACE  TEKPEBATUREt 
Development  of  Equations- 

In  calculating  the  maximum  surface  temperatures,  the  assump- 
liion  was  made  that  the  neutron  flux,  and  consequently  the  heat  flux,  was 
distributed,  as  in  the  KTR,  according  to  a  cosine  curve  which  dropped  to 
0,1^  times  its  maximum  at  the  end  and  the  plates.  Thus,  if  represents 

the  maximum  heat  flux  at  the  center  of  the  length  of  the  fuel  plate,  O.UQ^ 
represents  the  flux  at  the  end  df  the  fuel  plate.  In  addition,  it  was 
assumed  that  no  heat  was  conducted  linearly  along  the  fuel  plate. 

The  heat  flux  in  the  fuel  plate  can  now  be  expressed  as 

Q  =  Qq  sin  flx/L  (5) 

where  L  is  the  extrapolated  length  of  the  fuel  plate  that  has  zero  flux 
at  the  ends.  L  is  related  to  the  actual  length  of  the  plate,  a,  by 

L  -  a/  [l  -  (2)(0.1;ll6)/nJ  (6) 

-  a/0.7380 

Designating  the  temperature  of  water  as  T^  and  that  of  the  sur¬ 
face  as  Tj^,  the  water  flow  rate  in  pounds  per  foot  of  width  of  plate  per 
second  as  Gj  a  heat  balance  from  the  entrance  (designated  sub  1)  and  any 
point,  X,  along  the  plate  is 


3600  G  Cp  (T„  -  T  )  .  2  sin  trx/L  dx 


(7) 


Integrating  this  expression 


2  QqL 

T  «  T  +  " 

W  Wl  jm  TT 


(0.9165  -  cos  nx/L)  (8) 


Since  it  is  assumed  that  all  the  heat  generated  in  the  plate  is 
transferred  direetly  into  the  water,  the  following  relationship  is  valid  at 

any  point,  x. 

h(Tm  -  T„)  -  Qo  sin  nx/L  (9) 

Tjjj  »  +  (Qo/h)(sin  nx/L)  (10) 

Substituting  the  value  of  in  equation  (8)  into  equation  (lO), 
the  surface  temperature  may  be  expressed 


2  Q  L 


TfX  Q( 


T„  -  T,.„  +  ~  :~o: _ 

JSOOnCTG 


m  wi 


(  0.9165  -  cos  i^)  h"  L“ 


(11) 


In  order  to  find  the  point  at  which  the  maximum  value  of  Tj^  occurs 


'^'^m  ^  *^o  !E_  ^  Oo  «  5^  =  0  (12) 


or 

tan  TfX  ai  «  3600  Tf  Ct)  G 

L  2  h  L -  (13) 

The  value  of  the  heat  transfer  coefficient,  h,  was  obtained  from  an 
empirical  equation  given  in  McAdams.*  The  equation  is 


h  -5.6  (1  +  0.011t)(G')°'®/(D')°*^  (lU) 

Where  t  is- the  bulk  temperature  of  the  fluid  in  degrees  Fahrenheit,  O' 
is  the  pounds  of  fluid  flowing  per  second-square  foot  of  cross  section  and 
D'  is  the  hydraulic  diameter  in  inches.  The  diameter  was  taken  as  twice 
the  gap  width  which  is  approximately  equal  to  the  value  obtained  by  divi- 
ding  four  times  the  area  hy  the  perimeter. 


♦McAdams,  W.  H.,  o£.  cit.,  p.  183,  Equation  (9). 
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In  order  to  evaluate  this  heat  transfer  coefficient  and  also  the 
heat  capacity,  and  density,  it  is  necessary  to  know  the  average  bulk  tem¬ 
perature  of  the  liquido  This  is  taken  as  the  arithmetic  average  between  the 
inlet  and  exit  water  temperatures. 

The  total  heat,  H,  removed  from  a  foot  width  of  plate  in  one  second 
is 


(15) 


Integrating 

H  =  3.21;  X  10"^  Qq  i  Btu/ffc.*'s©c.  (16) 

The  temperature  rise  of  the  water  in  flowing  by  the  plate  is  found 

at  =  H/G  Cp  (17) 

As  a  first  estimate,  the  value  of  C  may  be  taken  as  1.00  Btu  per  pound- 
degree  Fahrenheit;  L  is  2,781  feet.  The  tenperature  rise,  in  terms  of  the 
maximum  heat  flux  and  flow  rate  is 

^T  =■  9.02  X  10"^  Q^G  (18) 


The  average  bulk  water  temperature  is  then  found 

T  +  (T  +  ^  T) 
t  “  1  1 


2 


(19) 


Method  of  Calculation- 

A  value  for  the  maximum  heat  flux,  ,  was  selected;  and  for  a 
series  of  values  for  the  mass  flow  rate,  G,  in  pounds  per  foot  of  width- 
second,  the  rise  in  water  temperature,  ^T,  was  calculated  from  equation 
(18).  Next,  an  entrance  water  temperature  was  selected  and  the  average 
bulk  water  temperature,  t,  was  calculated  from  equation  (19). 
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At  this  temperature)  the  values  for  the  density  was  obtained  as  the  recip¬ 
rocal  of  the  specific  volume  given  in  Keenan  and  Keyes#  and  the  heat  ca¬ 
pacity  was  obtained  from  Figure  12-8  in  OEML-156# 

For  the  known  gap  width  and  the  determined  density)  it  is  now 
possible  to  calculate  the  mass  flow  rate)  G',  and  the  linear  velocity,  V 
for  each  assumed  value  of  G.  Sufficient  values  are  now  known  to  calculate 
the  film  heat  transfer  coefficient,  h,  from  equation(lli). 

In  order  to  locate  the  point  of  maximum  surface  teii5)erature,  the 
values  of  tan  trx/L  are  calculated  from  equation(13))  and  the  corresponding 
values  of  sin  irx/L  and  cos  nx/L  are  read  from  tables.  Since  these  are  the 
only  values  needed  in  determining  the  maximum  surface  temperature,  it  is 
unnecessary  to  solve  for  the  quantity  x. 

The  values  obtained  are  substituted  in  equation  (11)  to  obtain  the 

maximum  surface  temperature. 

Ah  example  of  these  calculations  is  shown  in  Appendix  III, 

Results  of  the  Calculations- 

For  a  preliminary  group  of  calculations ,  an  average  bulk  water 
temperature  of  130^F  was  assumed  for  all  conditions,  and  a  density,  heat 
capacity,  and  heat  transfer  coefficient  were  determined,  based  upon  this 
assumption. 


♦Keenan,  J,  H,,  and  Keys,  F,-G,,  cit,,  p,  29-30 

♦♦Harrison,  W.  B,,  Physical  Properties  of  Water,  and  Water  Vapor^  ORNIi-l56, 
Part  12,  April  7,  19U9. 
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The  quantity,  maximum  surface  temperature  minus  inlet  water 
temperature  (T^  -  calculated  for  various  mass  flow  rates*  For 

each  of  the  proposed  variations  in  fuel  assemblies,  the  corresponding 
linear  velocity  is  calculated*  This  information  for  various  values  of 
maximum  heat  flux  is  given  in  Table  III  and  shown  graphically  in  Figure  2. 

More  exact  calculations  according  to  the  plan  outlined  in  the 
previous  section  were  made  for  the  19-plate  assembly  at  a  maximum  heat 
flvx  of  6  X  10^  Btu  per  square  foot-hour  and  inlet  water  temperatures  of 
100,  120,  11;0,  and  160^.  In  the  MTR  this  value  of  maximum  heat  flux 
corresponded  apf)roximately  to  a  total  power  of  30  megawatts*  The  results 
of  these  calculations  are  shown  in  Table  IV  and  are  presented  graphically 


in  Figure  3* 

A  calculation  was  made  to  determine  the  effect  of  changing  to 
an  18-plate  assembly,  maintaining  the  sai.xe  linear  velocity  between  the 
plates  and  increasing  the  heat  flux  proportionately*  The  overall  effect 
was  to  increase  the  maximum  surface  temperature  in  a  direct  relationship 
to  the  increased  flux*  This  amounted  to  from  3^  to  5^F  in  the  range 
of  velocities  from  20  to  1|0  feet  per  second,  $  to  7^F  in  the  range  l5  to  20 
feet  per  second  and  7®  to  lO^F  in  the  range  10  to  iS  feet  per  second* 

This  effect  is  shown  in  Figure  1;* 

In  calculating  these  values  of  maximum  surface  temperature,  the 
stipulation  is  made  that  boiling  does  not  occur  in  the  water  film.  Since 
this  is  the  desired  operating  condition  in  the  reactor,  these  results  are 
applicable*  If  the  flux  in  one  of  the  fuel  assemblies  is  distorted  by 
the  presence  of  a  shim  rod,  resulting  in  a  peaking  of  the  heat  flux,  a 
conservative  estimate  of  the  maximum  surface  temperature  is  found  merely 
by  substituting  this  maximum  heat  flux  in  the  derived  formulae* 


MAXIMUM  SURFACE  TEMPERATURE  MNUS  INLET  WATER  TEMPERATURE 
AS  A  FUNCTION  OF  WATER  FLOW  RATE  FOR 
VARIOUS  MAXIMUM  HEAT  FLUXES 
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MAX.  SURFACE  TEMP.— 


170 
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OPERATING  COiroiTIONS  FOR  VARIOUS  LIKITING  SURFACE  TEMPERATURES; 

Depending  upon  the  static  pressure  in  the  reactor j  various  niaxiitium 
surface  temperatures  may  be  allowed  without  incurring  the  risk  of  film 

For  any  temperatui‘e  the  minimum  allowable  water  velocity  corre¬ 
sponding  to  any  of  the  four  inlet  water  temperatures  can  be  read  from 
Figure  3  of  this  report.  I'or  these  flow  rates  the  pressure  drop  across 
the  fuel  assembly  can  be  read  from  Figure  2.  A  series  of  these  values 
are  given  in  Table  V  for  the  19-plate  assembly. 

RECO^affiNDED  CHANGES  IN  FUEL  ASSEMLIES? 

From  the  calculations  of  overall  pressure  loss  in  the  assembly,  it  is 
>  possible  to  isolate  the.  individual  pressure  losses  across  the  various  com¬ 
ponents.  This  study  illustrates  how  changes  might  be  made  to  reduce  the 

pumping  power,  ' 

For  the  19-plate  assembly  with  a  linear  velocity  of  30  feet  per  second 
through  the  plate  section,  the  exit  loss  from  the  lower  end  box  comprises 
approximately  of  the  total  pressure  loss  in  the  assembly.  By  changing 
this -exit  cross  section  to  a  2.2^0  by  2.2^0  inch  rectangular  opening,  this 
exit  loss  is  lowered  to  of  the  total.  Altering  the  upper  end  box  to 

a  rectangular  cross  section  would  restore  the  ability  to  invert  the  element 
and  produce  a  further  reduction  in  pressure  drop.  Eliminating  the  sudden 
change  from  circular  to  rectangular  cross  section  produces  an  additional 
6/?  reduction  in  pressure  loss. 


21t 


TABLE  V 

OPERATING  CONDITIONS  FOR 
VARIOUS  ALLOWABLE  SURFACE  TEMPERATURES  FOR 
19-PLATE  ASSEMBLY 


Allowable  Surface 
Temperature 
(°F) 

Inlet  Water 
Temperature 
(°F) 

Velocity  Between 

Fuel  Plates 
(ft/sec) 

Pressure  Drop 
Across  Fuel 
Assembly  (psi) 

210 

100 

23.7 

21.8 

120 

28.2 

30.lt 

lUO 

35.1 

1i6.3 

215 

22.3 

19,h 

120 

26.lt 

27.'0 

lUO 

32.5 

39.8 

220 

100 

21.1 

17.5 

120 

21i.7 

28.8 

lUO 

30.1 

3li.6 

160 

U0.2 

59.8 

225 

100 

20.1 

15.9 

120 

23.2 

21.3 

lUO 

28,0 

29.8 

160 

36.lt 

1^8.9 

230 

100 

19.0 

1U«U 

120 

21.8 

18.7 

liiO 

26.1 

26.3 

160 

33.1 

ltl.3 

235 

100 

120 

lUO 

160 

18.0 

20.6 

21t.5 

30.6 

13.0 

16.5 

23.lt 

36.8 
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A  sui!Tination  for  the  pressure  loss  in  the  assembly^  with  and 
without  these  inodifications  is  found  in  Table  VI 


TABLE  VI 

PRESSBEE  LOSS  IN  19-PLATE  ASSE^^BLY 
FOR  VELOCITY  OF  30  FT/SEC 
THROUGH  PLATES 


Pressure  Loss 

Pressure  Loss 

Without  Modifications 

With  Modifications 

Section 

C opponent 

(psi) 

(psi) 

Upper  End  Box 

Entrance  Loss 

3.86 

2.38 

Pipe  Friction 

0.1*8 

0.1*8 

E^^ansion  Friction 

0.11* 

0.11* 

Velocity  Head 

2.88 

2.88 

Expansion  from  Round 

0.81* 

— 

Fuel  Plates 

Entrance  Loss 

3.02 

3.02 

Friction 

12.61 

12.61 

Exit  Loss 

0.60 

0.‘'60 

Velocity  Head 

0,00 

0,00 

Lower  End  Box 

Sudden  Contraction 

1.39 

Friction  in  Pipe 

0.86 

0.86 

Exit  Loss 

10.92 

6.1*1* 

Velocity  Head 

2.88 

2.88 

Total 

3U.72 

26.53 
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APPENDIX  I 
NOMENCLA'I'URE 

A  Length  of  fuel  plate^  feet 
Cp  Specific  heat,  Btu/pound-^F 
Hydraulic  diameter,  inches 
D_  Equivalent  hydraulic  diameter,  feet 
f  Fanning  friction  factor 
g^  Gravity  constant,  32,2  feet/second 
G  Water  flow  rate,  pounds/foot  of  width-second 
G'  Water  flow  rate,  pounds/foot  -second 
h  Film  heat  transfer  coefficient,  Btu/foot^-hour-°F 
H  Total  heat  generated,  Btu/foot  of  width-second 

Kfs  Constant  relating  head  losses  due  to  sudden  contraction  to  flo^  rate 
K  Constant  relating  head  loss  due  to  expansion  to  flow  rate 
L  Lenght,  feet 
Ap  Pressure  loss,  pounds/inch^ 

Q  Heat  flux,  Btu/foot^-hour 
Qq  Maximum  heat  flux,  Btu/foot^-hour 
S  Specific  volume,  foot^ypound 
t  Average  bulk  water  temperature,  oF 
Tj^  Surface  temperature,  OF 

Temperature  change  in  water  flowing  past  fuel  plate,  ^F 
T^  Water  temperature,  °F 
u  Flow  rate,  feet/second 

Ujj  Flow  rate  in  smaller  cross  section,  feet/second 


Flow  rate  in  larger  cross  section,  feet/second 
V  Velocity  between  fuel  plates,  feet/second 
^  Density,  pounds/foot^ 

A  Viscosity,  pounds/foot-second 
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APPENDH  II 

CALCULATION  OF  PRESSURE  LOSS 

Calculation  of  Pressure  Loss  in  19-Plate  Assembly 
For  T  -  100°F 

/fe  0.72  C.  P.  »  0. 0001^814  /ft-sec 
A  62.0  /ft^ 

V  *  velocity  between  plates  in  ft/sec 
Flow  Areas? 

In  fuel  plate  section  . . . . . . 

In  2.250"  I.  D.  pipe  . . . . 

In  2.625”  le  D,  pipe  . . . . . 

In  3*019”  X  2.7U6"  area. . . . . 

Velocity; 

In  fuel  plate  section  . . . . . 

In  2.250"  I.  Do  pipe  . 

In  2,625”!.  D.  pipe  . . 

In  3.019  X  2.7U6»  area  . . . 

Above  and  below  assembly  . . 

Pressure  Loss 

In  upper  end  box- 

Entrance  loss, 

K  u^^ 

P  “  g  ' 


O.OliOli  ft^ 
0.0276  ft^ 

0.0376  ft^ 
0.0576  ft2 

V  ft/sec 
l.U^U  V  ft/sec 
1.07U  V  ft/sec 
0.701  V  ft/sec 
0.117  V  ft. sec 
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Friction  in  2,250"  I.  D,  pipe, 


T)  M  jO 

Ee  *  e  / 

A 


,  (2.25/l2)(l.U6U  V)(62.0)  =  3  ^2  ^  10^^  V 

(0.000U8U) 

2f  L 
‘  Ihk  D  gc 


=  2f  (1,4611  V)^(62.0)(7/l^ 

(2.250/l2)(32.2)(lUliy 


0.1783  V 


f 


V 

f 

p 

1 

0.0060 

0.00107 

5 

0.00U5 

0.0200 

10 

0,0036 

0,06ii2 

20 

0,0033 

0.235 

30 

0.0030 

0.1^81 

V  was  plotted  against  Ap  on  log-log  graph  paper,  and  the  resulting 
straight  line  was  found  to  be  of  the  form 
/IP  -  O.QOIOU 

Friction  in  5,50"  long  expander. 


AP 


2  gc  dUU) 


AP 


(0'.l5)(l.L6li  V  -  1.07k  V)^  (62,0) 
“  2(32r2)(lUU) 


0.000153 
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Expansion  from  round  section; 
_  \ 


(l.07)i  V  -  0.701  V)2  (62.0) 

2  (32:2)aiaD 


- 0.00093  V' 


Total  head  loss^  upper  box  - - —  0*00538  ^  +  0,0010^ 

Velocity  head  change 


^.79 


"  2  g„  dui; 


rl(0.701  v)^  -  (0.117  v)5  £62. gj -  0.00319 

^  202.2HlUUr 


p  170 

Net  pressure  loss,  upper  end  box -  0.0085?  V  +  0.00101^  V 

In  Fuel  Plate  Section- 


Entrance  loss. 


K 


AP 


2  duu) 


=  (0.g)(V)^(62.0)  _  0,0033U  V‘ 

2(32.2)(1UU) 
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Friction  in  plates, 


Re 


Du/® 


“  (0.23^)(V)(62.0) 

(12n0.000ii»l;) 


=  2.509  X  103  V 


2f  /O  L 

P  B  '  ' 

gc  ^  (II4U) 


_  2f  (62.0) (2.052) (12)  =  o  ftno  ^  v2 

'(32.'2)To:'23FKluiD 


V 

f 

.4p  ... 

1 

0.0100 

0.0280 

5 

0.0075 

0.525 

10 

0.006I; 

1.793 

20 

0.0055 

6.I6U 

30 

0.0051 

12.861 

1  79 

Again  this  relationship  is  fowad  to  bcj^p  *  - —  0,0291  V  * 

Exit  Loss, 

2  gcduu) 


Total  head  loss,  fuel  plates 


0.000598 


0.0039U  +  0.0291  7^*^^ 


Velocity  head 


0 


Net  pressure  loss,  fuel  plates 


0.0039U  V2  +  0.0291 
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In  lower  end  box- 

Sudden  contractionj 


K 

2  g^dhk) 


(0.20) (0.701  V)^  (62.0) 
2(32.2)  dhU) 


0.00065? 


Neglect  grac^ual  contraction* 

Friction  Drop,  taking  same  Re  as  in  upper  end  box. 


2f  u^y^L 

gc  dUuu) 


-  2f  (I.ii6u)^(62.0)(l2.5)(l2) 
(2.25)(32.2)UU1+K12) 


0.318  f 


V  1 

f 

-  p . 

1 

0.0060 

0.00191 

5 

0.00U5 

0.0357 

10 

0.0036 

0.111/? 

20 

0.0033 

0.ii20 

30 

0.0030 

0.859 

This  is  found  to  beZlp  «  -  0.00186  V 


Exit  loss, 

.  ‘"H  - 

=  (l.li6li  V  -  0.117  V)^(62.0) 
2(32.2)(iltU) 

Total  head  \loss,  lower  end  box  - 

Velocity  head  (negative  of  upper) 
Net  pressure  loss,  lower  end  box  - 


2 


-  0.01213 

•—  0.01279  v2  +  0.00186 

- - - -  0.00319 

0.0096  +  0.00186 
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Suimnaiy  of  pressure  losses: 


In  upper  end  box  - — —  0*00857  +  O.OOlOIj. 

In  fuel  plates  -  0.0039U  +  0,0291 

In  lower  end  box  - — —  O.OO96O  +  0.00186  V 


Total 


0.02211  + -0.0320 
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APPENDIX  II 

CALCULATION  OF  KUXIMUK  SURFACE  TEMPERATURE 

For  a  maximum  heat  flux  of  6  x  10^  Btu/ft^-hr,  the  total  heat  transferred 
from  a  foot-wide  strip  of  fuel  plate  is 

H  =  3.2iil  X  10"^  Qo  L 

=  (3.2U1  X  10"il)(6  X  10^)(2.781) 

V 

-  %0.1  Btu/ft-sec 

For  a  heat  capacity  of  1.00  Etu/j^°F,  the  rise  in  the  temperature  of  water 
in  flowing  past  the  fuel  plate  is 

m  -  H/G.Cp 

=  5U0.7/G 

G  is  the  mass  flow  rate  of  water  past  a  foot-wide  strip  of  plate.  This 
value  may  be  related  to  the  mass  flow,  G',  per  square  foot  of  flow  area  by: 
for  gap  width  of  0.1177*  =  0.009808' 

G'  =  G/O. 009808  =  102.0  G 
The  average  bulk  water  temperature,  t,  is  found 

^W]_  +  +4t) 

t  - - - - 

2 

The  inlet  water  temperature  is  100°F.  The  film  heat  transfer  coefficient,  h 
is  found  by 

‘  h  »  5.6  (1  +  0.011  t)(G')°*®/D°*^ 

The  following  quantities  may  now  be  tabulated: 


35 


G 

(?j^ft-sec 

AT 

(°F) 

G 

0)*/ft^-sec) 

-PO 

h 

(Btu/hr-ft2-°F) 

6 

90.1 

612.0 

11;5 

3788 

8 

67.6 

816.0 

13U 

U531 

10 

%.l 

1020,0 

127 

5261 

12 

U5.1 

1221;.  0 

123 

5961 

lU 

38.6 

1|;28. 

119 

6633 

16 

33.8 

1632. 

117 

7301 

18 

30.0 

1836 

115 

79U8 

20 

27.0 

201;0 

i;)4 

8566 

22 

2i;.6 

1 

221J; 

112 

9206 

The  point  of 

maximum  surface  temperature  is  noi 

flT  investigated. 

The  values  ^ 

are  determined  for 


tan  Hi  .  2033  G 

1*  "  “h — 


G 

Ol^^ft-sec) 

Tan  M 

Xi 

Sin 

L 

^  TTX 

Cos  L 

6 

-3.2202 

0.9550 

8 

-3.5895 

0.9633 

-0.2681; 

10 

-3.861;3 

-0.2505 

12 

-U.0933 

O.97IU 

-0.2373 

114 

-li.2907 

0.9739 

-0.2270 

16 

-i;.);556 

0.9757 

-0.2190 

18 

-1;.6037 

0.9772 

-0.2123 

20 

-1;.71;67 

0.9785 

-0.2061 

22 

-1;.8588 

0.9795 

-0.2016 

Now  the  maximuin  suri’ace  temperature  may  be 


found  from 


■m 


T  + 


2Qo  L 


p « 


(0.9165  -  cos  22L)  +  5o  Hi 

L  h  L 


=  100  +  6  X  10 


5 


r (0.9165) 

L  20337 


cos  nx 


h  G 


+  1  sin  M 

h  L 


3<5 


The  inner  velocity  can  be  found  from  the  mass  flow  rate  and  density 
(reciprocal  of  specific  volumes) 


V 


»  ^  102.0  G  S 


These  are  found  tc^  be 


G 

(#/f  t-sec  ) 

T 

m 

(Op) 

H 

6 

311 

0.01632 

10*0 

8 

271 

0.01626 

13.3 

10 

2U3 

0.01623 

16.6 

12 

226 

0.01622 

19.9 

lu 

212 

0.01620 

23.1 

16 

201 

0.01619 

26. 

18 

192 

0.01618 

29.7 

20 

185 

0.01618 

33.0 

22 

179 

O.OI6I8 

36.3 
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